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Recent advances in nanotechnology and atomic physics may allow for a demonstration of the 
dynamical Casimir effect. An array of film bulk acoustic resonators (FBARs) coherently driven 
at twice the resonant frequency of a high quality electromagnetic cavity can generate a stationary 
state of Casimir photons. These are detected using an alkali atom beam prepared in an inverted 
population of hyperfine states, with an induced superradiant burst producing a detectable radio- 
frequency signal. We describe here the results of the simulations of the dynamics of superradiance 
and superfluorescence, with the aim to optimize the parameters for the detectability of Casimir 
photons. When the superradiant lifetime is shorter than the dissipation time, we find superradiant 
evolution to be similar in character but dramatically slower than in the usual lossy case. 

PACS numbers: 12.20.Fv, 42.50.Pq, 85.85.+j, 42.50.Lc 



I. INTRODUCTION 



Observable effects due to the change in the boundary conditions of quantum fields, like the creation of particles in 
an expanding universe [l[ or the Casimir force Q , provide crucial information on quantum vacuum at the macroscopic 
level. After the recent results on Casimir forces, with measurements performed in a variety of geometries ranging 
from the original parallel plane [1, H[ to the sphere-plane [f| @, 0, H, [1] and crossed-cylinders 
to understand dissipative effects of vacuum fluctuations, especially its interplay with relativity 
dissipation mechanism should induce irradiation of photons, a phenomenon also known as dynamical Casimir effect 
EE, EE!]- This can be understood both as the creation of particles under non-adiabatic changes in the boundary 
conditions of quantum fields, or as classical parametric amplification with the zero point energy of a vacuum field mode 
as input state. In this paper, following on the proposal described in [l9j|, we describe a model for the superradiant 
amplification scheme with particular emphasis on its dynamics and the optimization of the involved parameters. 
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II. GENERATION AND DETECTION OF CASIMIR PHOTONS 



As discussed in more detail in (l9| . under parametric amplification in an electromagnetic cavity an initial state of 
Nq photons with frequency within the resonance bandwidth of the fundamental mode of the cavity uj is transformed 
into a squeezed state with an average number of photons growing in time as 0, [l?], EH] : 

N C3s (t) = N sinh 2 (Lj mcch et), (1) 

provided that the parametric resonance condition with a mechanical driving at a frequency w mec h = 2o> is fulfilled. The 
term u) mec hC in the hyperbolic sine function represents the squeezing parameter, with the modulation depth e = v/c, 
where v is the velocity of the resonator and c the speed of light. This exponential growth is eventually limited by 
the photon leakage of the cavity expressed through its quality factor Q, which saturates at the hold time r = Q/to, 
reaching a maximum photon population: 

A^cT = ^cas(r) - A^o sinh 2 (2Qe). (2) 

Given an initial number of photons in the cavity, which will be attributable to the quantum vacuum at temperatures 
such that KbT <C fno (see also [2(| for a detailed discussion), the average number of photons at saturation in Eq. @ 
strongly depends on the product of two parameters, Q and e, which can be on the order of 10 8 and 10 -8 respectively. 
The expected saturated power initiated by Casimir emission is: 

Pes = N%?— (3) 

T 

and for a 3.0 GHz FBAR resonator and Qe ~ 1, the saturated power is 3 x 1CP 22 W, which is too low to be detectable 
using current technology. This demands the use of an efficient, nearly quantum-limited, photon detector in the 
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radio-frequency range. Ultra-sensitive atomic detection schemes can be exploited for detecting Casimir photons by 
preparing an ensemble of population-inverted atoms in a particular hyperfme state, which for alkali atoms ranges 
from 0.2 GHz for Li to 9 GHz for Cs, whose transition frequency corresponds to the cavity resonance. An additional 
amplification process is available in which the weak Casimir signal triggers the stimulated emission of the ensemble 
of atoms. This effect is a form of superradiance [2ll. I22I] . The hyperfme transition in the ground state occurs through 
a magnetic dipole interaction, and its natural lifetime in free space is approximately: 

^ 4tt 3ft 

where \ib is the Bohr magneton and /xo the magnetic permeability in vacuum. This natural lifetime in free space is 
favorably reduced inside a resonant cavity due to the modification of density of states [HI, [24[ : 

T cav _ 47T 2 V 47T HV 

1 3Q A 3 ^mqSttmIQ' [ ' 

where V is the cavity mode volume. For a few GHz cavity with 1 cm 2 cross-sectional area and Q = 10 8 , the natural 
lifetime is reduced by a factor of 10 10 . In spite of this cavity-enhanced spontaneous rate, the typical hyperfme 
transition lifetime for the alkali atoms is still impractically long, on the order of 10 3 -f- f 5 s. To shorten this timcscalc, 
let us suppose to inject A at excited atoms into the cavity. The Casimir field acts on all atoms, stimulating emission on 
a time scale of the supcrradiant lifetime, defined as Tsr = Tf av /iV a t, which is in the millisecond range for A at « f 8 
or less. An atomic density large enough will then induce a supcrradiant burst with peak power of Psr ~ N a thu>/Tsn, 
increasing quadratically with the number of atoms. Considering a few GHz resonator with 10 8 atoms and Tsr = f 0~ 3 
s, yields Psr = 10 -13 W, a billionfold improvement over the power without superradiant amplification as in Eq. ([3]). 

Spontaneous emission into the cavity mode by the atoms will also trigger a superradiant burst, a process also known 
as superfluorescence. To distinguish this source of background from the Casimir stimulated superradiance signal, one 
may study the temporal intensity envelope of the amplified photons. Both the average delay (Td) of the peak intensity 
from the initial excitation of the atomic population and its fluctuation (ATd) are reduced with increasing number of 
resonant photons -/V p h initially present [25[: 

N K 



T D = Tsr In * , AT D = 2T SR / y/1 + N ph . (6) 
\ i + iVph / 

It should be noted that these standard results follow when the cavity lifetime is much shorter than the supcrradiant 
lifetime. For very high Q, the system of coupled equations of motion must be integrated directly. Tailoring the atomic 
number can further distinguish the Casimir stimulated superradiance from superfluorescent pulses. In order for the 
superradiant pulse to develop fully, the growth rate must exceed any decay process, which is primarily due to Doppler 
dephasing in the atomic cloud, and the atoms must remain in the interaction region for a time longer than the delay 
time. A proper choice of A at may suppress superfluorescence relative to Casimir superradiance provided that the 
atoms will be removed from the cavity after the expected Casimir delay time but prior to the superfluorescence delay 
(■Wph = in Eq. i[6])). The supcrradiant emission can be detected by coupling a power or field detector to the cavity, 
with response time small enough to resolve one superradiant lifetime. One issue with this direct measurement is the 
possible reduction of the quality factor of the cavity especially for large coupling efficiency. Micro-bolometers mounted 
on etched "spider-webs" have an ultimate sensitivity of 10 _16 W/\/Hz in the GHz range [26[. Spectrum analyzers 
are sensitive to sub-fW RF power of kHz bandwidth (2?| , and the temporal profile of the burst can be reconstructed 
through vector analysis. 



III. SUPERRADIANT AMPLIFICATION MODEL 



To validate these qualitative estimates, it is necessary to simulate the dynamics of an atomic beam travelling 
through a low dissipation cavity. We shall discuss here specifically Na atoms initially optically pumped into the 
\F = 2,mp = 2 > ground state, though the approach is general. Let us consider the second-quantized Hamiltonian 
with relativistic and hyperfme corrections absorbed into the unperturbed atom term H atom [28| : 

N 2 

H = Y, H T° m - ZT.VJ ■ Mr 3 ) + |A( rj )| 2 - to ■ Bfo) + £ f^ x (a{a x + 1/2) . (7) 

j=l ] J A 
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FIG. 1: Detectability plot in the e — Q plane. Depicted from top to bottom are the regions in which direct detection is possible 
with state-of-art radiofrequency detectors, the region where superradiant amplification extends the range of detection, and the 
region experimentally inacessible also to superradiant amplification, with its lower limit due to the speed of micro-bolometers 
and heterodyne receivers of the current generation. The dashed line corresponds to the benchmark values of Q and e satisfying 
Qt = 1 which are at the edge of current technology of superconducting cavities and nanotechnology. 



The subscripts j and A are the atom and field mode index, respectively, lo\ is the mode angular frequency. The fields 
A(r\j) and B(rj) = V x A(iVj) are defined in terms of field profiles functions U(rj) and creation and annihilation 
operators (a,aJ) as: 



A(r,) = £ 



2e ui\ 



Bfo) = VjM)ftu>j2 k(ij) x U(r i )a A - k(r?) x U(r,- 



(8) 
(9) 



where the field profile functions form an orthonormal set and the a,a^ operators fulfil the usual commutation rela- 
tionships: 



d 3 rU A (r,)U x ,( ri ) 



0. 



(10) 



Given no initial population on an upper state of an electric dipole transition and no initial or applied field resonant 
with the same, the second term can be ignored and each atom evolves only within the manifold of ground hypcrfinc 
states. Also, the third term is negligible compared to the fourth in this case. The atomic Hamiltonian then can 
be represented by an 8 x 8 matrix and the Hcisenberg equations of motion derived with some effort. To simplify 
the discussion, let us assume a cylindrically symmetric cavity with field propagation (k) primarily along the cavity 
axis (z), with quantization axis along the cavity axis, so that the active modes will be circularly polarized. If 
we also assume the ideal case where the Na atoms are prepared in the \F = 2, nip = 2) state, then only the 
\F = 2, mp = 2) — \F = 1, mp = 1) transition will be active and the Hamiltonian can be reduced to: 



3=1 







using the standard notation for polarization, 

\x = /!+£+ 



) + huJx ( a i 



\ax + 1/2 



[is- + (i z z, 



(11) 



(12) 
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and likewise for B and U with e± = (x±y)/\/2. The hyperfinc transition resonance frequency ft = {H^f" 1 — Hfl° m — 
Hz =2 B® (vj) + B® c (rj)} /H can be tuned with an applied DC magnetic field along the quantization axis. Fur- 
thermore, given = /i_ = (U, it is convenient to define the Rabi frequency coefficient X\( r j) = ^1/ pbofkJ\/2 U\{rj)/h 
and employ Pauli matrix operators in writing the Hamiltonian as: 



H = ^ hw x [a{ a x + 1/2J + (13) 

A 



3=1 k A 



Here it is understood that the field operators represent only right circular polarization and the atom operators act 
only on operators of the same atom. The Hcisenberg equations of motion resolve to: 



{xA(rj)o- +i 6 A +XA( r i)^-34} ( 14 ) 

A 

a+j = i(n-k- Vj)o-+j +X]^( r 3)^3«A ( 15 ) 

A 

a A = -i^AflA + ^^A( r j)o'-i- (16) 

3 = 1 

Note that total energy, excitation number are conserved in a lossless system as well as the Bloch vector length. Optical 
losses from the cavity can be included by adding a negative imaginary loss rate, T = uj\/2Q, to the mode frequency. 
These equations suggest the operators can be factored as a\ — a\e~ lUxt and a±j = cr±je =t *( f2 ~ k ' v ^ t , isolating the 
variable envelopes which are easier to integrate numerically. Furthermore, if the cavity lifetime is much shorter than 
the superradiant lifetime (1/T <C Tsr), then it is convenient to first integrate Eq. (fH)|) approximately through a 
Green's function method by assuming the atom does not evolve significantly within one cavity lifetime [281 ]: 

N 

axMnaxWe-^+^ + J^xiir^gxiWZ-j, for rT Sfl »l, (17) 

3=1 

where g\j{t) = (1 — exp[(— iA\j — T)t\) / (iA\j + T) with detuning between cavity and atom resonance is defined as 
A\j = lo\ — (fi — k • Vj). Essentially, each atom responds to the emission from all atoms accumulated over one cavity 
lifetime. This result can then be inserted into Eqs. p4[) and p5[) . defining the detuning between atom resonances 
Aij = k • (Vj - Vj), to yield: 

~a zo = -2^ jxA^-^^W (18) 
+ ^{^^^A^O^We-^^+^-. + XA^OxA^OsA.We^" 4 ^^} \, 



^+3 = E xl(rM~ lAx >- V)t °zMQ) +J2^ r M r ^xS)e iAi ' t a zj a. i . (19) 

A I t=l ) 

The field intensity, square of Eq. (fTF)) . can then be produced from the atomic state by first numerically integrating 
these reduced equations. The terms linear in a\(0) represent the negligible action by the initial field. The quadratic 
terms involve feedback from all other atoms that generates superradiance. Note that on resonance (A\j = 0) we 
arrive at the expected time scale, for t <C 1/T, through the coefficient: 

xKrMrMt) = ^f^^A UliTj) U x {v t ) = Ufa) U x (r % ). (20) 

For significant amplification, the number of interacting atoms must be large and it becomes prohibitively difficult 
to solve the 27V coupled equations of motion. Instead one can dice the atomic distribution in phase space into finite 
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FIG. 2: Dynamics of superradiant amplification of Casimir photons (assumed in the number of 10 2 ) and competition with 
superfiuorescence. On the left plot, the superradiant emitted photons (solid line, black) compared to the superfluorescence 
in the absence of Casimir photons (dotted line, red) are depicted versus time. The dot-dashed curve in green indicates the 
field strength profile felt by the atoms. On the right plot, the atomic inversion populations for superradiant (solid line, black) 
and superfluorescence (dotted line, red) are also plotted versus time. The dot-dashed curve curve in green indicates the 
degree of conservation of the Bloch vector length during the simulation. The simulation is performed for for Na atoms with 
Q = 10 9 ,iV ot = 10 10 , 1 m/s atom velocity, 10 mK temperature. The relevant time scales are: 1/T = 0.18s, Tf av = 10 6 s, 
Tsr = 10 _4 s, and Doppler dephasing time TJ = 63 ms. 



close-packed cells so that the number of cells is computationally feasible and the pertinent distribution and field profile 
features are retained. These constraints typically imply that the number of atoms per cell is at least ten. In this case, 
the correlation between the atomic population and photon number is inconsequential (i.e. atom stimulation rate is 
not significantly limited by the stored energy) and the expectation value of the operator products can be factored as 

< a+j&x >~< d+j >< a\ >. The sum over atoms can then be replaced with the multiple sum over cell indices: 

TV 

E /i( r i. v i) -> E E E E N w m< x < y < z < v \\ >? 011 )' ( 21 ) 

j=l i j k l 

provided the atom position and velocity are replaced with the mean cell position tj — >< r > ce u and velocity Vj — >< 
v > ce u and the atom operators now represent the typical atom in the cell. The number of atoms in cell {i,j,k,l} 
is designated as Nijki, and only the parallel component of velocity is relevant. The cellular equations of motion 
can then be numerically integrated given the initial expectation values, atom and field distribution, atom resonance 
and moment, and field resonance. Collisions can be simulated by introducing a random phase factor, depending on 
the number of atoms in the cell, to the atom operators between time steps. Assuming an initially inverted atomic 
population, superfluorescence can be simulated by assuming zero photons present and an initial "tipping angle", 

< a Z j > (t = 0) = 2/y/Nat, of the Bloch vector representing spontaneous atom state fluctuations [2a]. Likewise, 
Casimir stimulated superradiance may be simulated by assuming the same tipping angle and initial field amplitude 

< a\ > (t = 0) = y/N™^, given the parametric amplification process generates a coherent state of the field. 



IV. DISCUSSION OF NUMERICAL RESULTS AND CONCLUSIONS 



The coupled equations of motion, Eqs. (|14 [) -(|16 [) . were integrated by Rungc-Kutta method for a small sample of 
Na atoms traversing a planar-concave optical cavity through the diameter at 10% of the cavity length height above 
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TABLE I: Configurations yielding effective discrimination between Casimir superradiant signals and superfluorescence noise. 
The degree of discrimination is indicated by the ratio of the peak intensity for the SR case to the SF case (r\) and the analogous 
ratio of the ground state population after the atoms exit the cavity (£). The ground state population of the atoms near 
resonance is given by p^ d = (1 — (o"z(A ~ 0)))/2. The superradiant lifetime for these cases is in the range of Tsr = 1 — lOO^is, 
roughly 10 3 times less than the cavity decay time, implying FTsr <C 1. 



the planar mirror. The cavity is tuned so that the first longitudinal TEMOO Gaussian mode is resonant with the 
atoms. The atom distribution was divided into five cells varying in velocity symmetrically about the mean speed. 
The number of atoms, cavity quality, atom mean speed, atom transverse temperature (i.e. velocity spread), and initial 
photon population were varied to identify promising configurations and sensitivity of the superradiant output. We 
found that there is a range of configurations, corresponding to the high end of experimental feasibility, that achieve 
a high degree of discrimination between the Casimir stimulated superradiant and the superfluorescent emission. In 
Fig. 2 we show the typical emission profile (left panel) and atomic inversion population (right panel). As suggested 
above, the superfluorescent emission is suppressed if the atoms leave the field mode region before its amplification 
can develop. Even a single photon in the cavity can generate a signal 60 % above the superfluorescent background. 
Atom numbers too high or too low, depending primarily on Q and atom speed, allow both cases to develop fully 
or very little, respectively. Superradiant amplification is quite effective for sufficiently high gain, casing RF detector 
sensitivity requirements. The gain is very sensitive to N at and Q and the atom transit time (oc l/w at ). Also, the 
strong discrimination in amplification indicated by r\ in Tabic 1 implies that the detector response time need not 
be a significant constraint, as would be the case if the delay time should be measured. The discrimination ratio is 
fairly insensitive to temperature, suggesting that a slowed beam rather than a trap could be used, providing larger 
N a t- Strong amplification necessarily results in a large fraction of the atoms driven to their ground state. Therefore, 
probing the atomic beam ground state through D-line excitation after it transits the cavity, as we had proposed in [il^ ] , 
is a viable option that avoids detector losses in the cavity. The high values of N at and Q would be eased significantly if 
the number of initial Casimir photons is increased. Equation ^ implies N^^. = 0.7, 6, 370, 10 6 for Qe = 0.5, 1, 2, 3, 4. 
The rapid change in initial photon population due to Q for a given modulation mechanism is not taken into 

account in Table 1. One caveat to the above results and discussion is that the computed time for the superradiant 
burst to develop is far longer than the standard prediction, Eq. ([B]). The configurations tested imply TT$r -C 1, in 
violation of the approximation condition under which Eq. ^ was derived. Indeed, the reduced equations of motion, 
Eqs. (fl~7|) and (fT9]l . were quite unstable under these conditions. As we have shown, the standard result develops 
naturally from the general theory given a lossy environment. The transition to a low loss environment appears to 
have a dramatic effect delaying the field growth. We arc currently studying the origin and domain of validity of this 
effect. 
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